Manipulating single enzymes by an external harmonic force 
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We study a Michaelis-Menten reaction for a single two-state enzyme molecule, wliose transition 
rates between the two conformations are modulated by an harmonically oscillating external force. 
In particular, we obtain a range of optimal driving frequencies for changing the conformation of 
the enzyme thereby controlling the enzymatic activity (i.e. product formation). This analysis 
demonstrates that it is, in principle, possible to obtain information about particular rates within 
the kinetic scheme. 



Recent advances in single molecule spectroscopy have 
made it possible to follow the catalytic activity of indi- 
vidual enzymes over extended periods of time [J, [3, S, 
0, H, 0| • Although the catalytic rates have been shown 
to be distributed in time, the general Michaelis-Menten 
scheme, originally proposed for an ensemble of enzymes, 
turned out to hold quite well also on the level of a sin- 
gle enzyme Q, i, i, i, i, @, 0, Hi- What is still left 
open is the nature of the conformational changes and 
its relationship to the enzymatic activity @, 0, S l^- 
Deeper insights into the conformational-activity relation- 
ship can be obtained by influencing single enzyme con- 
formations by an external force and following the effects 
on the catalytic activity. The possibility to manipulate 
enzymatic turnovers by applyinga mechanical force has 
been recently demonstrated in [105. Also, voltage-gated 
ion channels can be manipulated electrically to switch 
between two states ll|. Here, assuming a two-state 



Michaelis-Menten scheme for a single enzyme, we inves- 
tigate how one can manipulate the activity (turnovers) 
of an enzyme by an external harmonic force. The oscilla- 
tory force is shown to be able to increase the weight (oc- 
currence time) of a desired conformation in the scheme, 
controlling the enzyme activity. 

The reaction scheme that we examine is shown in 
Fig. [H We model the effect of the time-dependent ex- 
ternal force as a shift in the conformational energy land- 
scape. Assuming that the transitions between different 
conformations, with rates a^- and (3ij, follow an Arrhe- 
nius activation, and using a harmonic force which mod- 
ulates barriers separating conformational states, we find 
the following time dependence of the rates 
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Here, lo is the angular frequency of the external driving 
force, the initial phase of the force; and q;|°\ (3^^\ 

e^"^ , and e\^'^ are phenomenological constants character- 
izing the enzyme and the amplitude of the driving force. 
In Eqs. dl]) we assume that the potential energy land- 
scape of the conformations responds instantly to changes 
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FIG. 1: Michaelis-Menten scheme for a single enzyme with 
two conformations. E denotes the enzyme molecule, S the 
substrate, and P the product. Rates aij and j3ij {i,j = 1, 2) 
stand for transitions between conformational states 1 and 2, 
while rates k^i {x = —1, 1, 2) quantify internal conversions. 



in the force. This assumption can be relaxed, for instance 
by including different phases (/jq for aij and Note, 
however, that a 180° phase difference, corresponding to 
a change of sign for e-^\ is included in Eqs. H]). We as- 
sume that the substrate S is in great excess, such that the 
rates for the reaction with the substrate like ku = fcj°' [S] 
become independent of time. The scheme in Fig. [1] is 
closely related to schemes describing resonant activation 
in discrete systems fl^ . 

Simulations. As shown below, for certain limits of the 
rates analytical expressions can be derived. To obtain in- 
sight into the general behavior of enzymatic activity in- 
cluding noise, we use a stochastic simulation. To this end, 
we implemented the Gillespie algorithm that pro- 
vides random values for the waiting time between two re- 
action steps, and the direction of the reaction as weighted 
by the corresponding Arrhenius factor. We start the sys- 
tem in one of the states E^ -I- S, and then perform jumps 
to one of the neighboring states according to the Gille- 
spie reaction probability. Which state it jumps to is de- 
cided by picking a waiting time for jumping to state ESi 
according to the exponential distribution with rate pa- 
rameter ku, and one for jumping to Ej -I- S according to 
the cumulative distribution 
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and then the jump is made to the state with the shortest 
waiting time. Subsequent jumps to neighboring states 
are determined similarly, and one iteration of the sim- 
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FIG. 2: Mean turnover time (r) as function of driving fre- 
quency ijj for the aij,l3ij. Horizontal lines repre- 
sent the limit cases derived in the text, (i) ku = k-u = 1, 
k2i = 0.3, k22 = 3, a\l^ = 80, a$,°' = 800, p'^^ = 160, 
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lowing parameters were changed: ku — fc_ii = k-12 — 0.01, 
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£21 = 2 (units on right axis). Simulations 



were run for 5 x 10 turnovers per frequency. 



ulation stops at time t = t when the system reaches 
one of the product states E.^ -I- P. The conformation i 
of the enzyme is then used as the initial state for the 
next iteration with the initial phase i^q'^''^'' = ujt + 0q°''^'' 
determining the new momentary force. Performing itera- 
tions over many periods of the oscillating force we obtain 
a time series for the turnover times r, from which the 
average waiting time (r) yields. Plots of (t) as a func- 
tion of oj are shown in Figs. [5] and [51 Depending on the 
choice of parameters, switching between two or three ki- 
netic regimes of the system is revealed. Certain limiting 
behaviors are accessible analytically, as discussed soon. 
For the parameters used throughout this work the rate of 
product formation in the conformation 2 is higher than 
in conformation 1. 

We note that in our analysis the parameters (except for 
case (ii) in Fig. [J) are chosen such that detailed balance 
is fulfilled, as required for certain single enzymes that are 
coupled to the surrounding heat bath 1J|. The explicit 
condition arising when detailed balance is applied around 
the loop in the reversible reaction pathway in Fig. [1] is 
k-iiP2iit)ki2ai2{t) = kiiPi2{t)k^i2a2i{t). The condi- 
tion is applied for any momentary value of the external 
force, and therefore any time t, since the system could 
equally well be held constantly at these forces. Detailed 
balance could be violated, for instance, by enzymes con- 
verting ATP energy during their cycle. 

Fast switching limit. Consider the case when the 
switching between enzyme conformations is much faster 
than the reactions with the substrate, as displayed in 
Fig.H Note that case (ii) in Fig. [2] has different fcn's [Hi. 
It also includes a 180° shift in the phases of variation of 
'ij under the external force distinguish- 



FIG. 3: Plot of (/(<^)) as a function of (p for different driving 
frequencies lo, for the same conditions as in Fig. [5] (i). 



ing between substrate-bound and substrate-free enzyme 
states. Parameters were chosen to reflect situations with 
enhanced enzyme efficiency at intermediate frequencies. 
Different scenarios, for instance with decreased efficiency, 
are also possible. In both cases in Fig. [2] the rates satisfy 
kxi ^ ctij , (3ij . This means that we can define eflFective 
rates k^ as averages over the various conformations. The 
weights in these averages are given by the probabilities 
for the enzyme to be in the different conformations. If uj 
is much faster or slower than the relaxation time between 
the two enzyme conformations, these probabilities can be 
easily obtained. For instance, the probability for the en- 
zyme to be in conformation Ei, given that no substrate 
is bound, is 
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where (j> € (0, 27r) is the phase of the force, i.e., ujt + (j)o = 
<f> + 27rn for some integer n, and a bar over a quantity 
means that this quantity is averaged over all phases. 
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where Iq is the modified Bessel function of order 0. The 
effective rates kx can now be obtained in the two limits. 
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where x in p{i\x,(f)) refers to the unbound state E + S 
for X — 1, and to the bound state ES for a; = — 1 and 
X — 2. The rate of the overall reaction is then given by 
the standard Michaelis-Menten expression [3, 0, 13| 
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as function of the phase (j). The average waiting time 
(t) for one turnover can now be found by averaging over 
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FIG. 4: Estimate of p{2\(ji) as a function of for the same pa- 
rameters as in Fig. |5] (i) . Note that there are only two different 
analytic limits, since two of the three limits are identical for 
the chosen parameters. 



FIG. 5: Plot of the fraction p{2) the enzyme spends in con- 
formation 2 (lower curve), as a function of driving frequency 
uj for the same parameters used in Fig.[2](i). The upper curve 
depicts the fraction of product formation while the enzyme is 
in conformation 2. 



all phases and calculate (r) = The values of (r) 

calculated according to Eqs. (O and (O are shown in 
Fig. [2] by the horizontal lines. We note that the limits 
where u) <^ a^- , /3ij can also be obtained as limit cases of 
the system studied in Ref. [3]. 

Another experimentally relevant quantity is the tem- 
poral probability density for the formation of products 
at a given value of the force as determined by the phase 
(f). A plot of this intensity of product formation events, 
which we label by is shown in Fig. [3l An ap- 

proximate bimodal behavior is found for certain choices 
of the parameters. To obtain the analytic limits we 
note that if the external force is slowly varying, i.e., 
^ ^ kxi,ctij,Pij, then the intensity will simply be equal 
to the rate (/(</>)) = (</>)■ This expression also holds 
for Lo ^ aij , Pij where it is independent of (j). For the 
intermediate cu where k^i ^ cj <C aij , Pij it will be the 
distribution of conformations in the final step of the re- 
action and the corresponding rate constants that deter- 
mines when the product formation happens, such that 
{I{4>)) = X^i ^2iP(* |ES, (?!))/fc2- We point out here that 
at slow frequencies it is the overall efficiency of the en- 
zyme that is probed at different magnitudes of the ex- 
ternal force, while at the intermediate frequencies it is 
only the relative efficiency of the final step in the reac- 
tion which is probed at different external driving. Thus 
tuning the frequency allows one to selectively study the 
final reaction step. 

The behavior of {I{<f))) in Fig. [3]can be compared with 
the probability of being in, say, conformation 2 (the fast 
one) as a function of (/>, as shown in Fig. 01 Note the 
large change in amplitude for the chosen parameters. To 
obtain the analytic limits we write the probability as 
p{i\(j>) = p(i|E,0)p(E|0) -hp(«|ES,(A)[l -p(E|0)], where 
p(i|E, (j)) is given by Eq. ([3]). To frnd p{E\(j)) note that at a 
given phase (f> the effective rate for the system to leave the 
state without bound substrate is fci(0), and the rate for 



leaving the state with bound substrate is k-i{(/)) + k2{4>)- 
In fact the quasi-steady-state probability for being in the 
state without bound substrate is 



fc_i(</))-ffc2W 
fci(0)+fc_i(0)+fc2(<^) 
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in the two cases given by Eq. ([5]) where either lu <^ k^i or 
S> kxi- The resemblance of Figs. |3] and [5] reflects the 
dominating role of the conformation 2, due to its larger 
catalytic rate, in the enzymatic activity. 

To illustrate further the mechanisms behind the tran- 
sitions of (r) with changing driving frequency we plotted 
the fraction of time the enzyme spends in conformation 
2, labeled p{2), during a simulation run in Fig.O The an- 
alytic limits of p(2) can be obtained straightforwardly by 
averaging the probability of being in conformation 2 over 
all phases, p{i) — J^^ [dtj) / {2t:)] p{i\(t)). This behavior 
can be compared with the fraction of product formation 
events that occurs while the enzyme is in conformation 2, 
labeled by p{2\E + P) in Fig. [5l The analytical hmits of 
p{2\E + P) can be obtained by averaging over all phases 
<f) the effective rate constant for forming a product in con- 
formation i, fc2ip(i|ES, (j)), divided by the overall effective 
fc2((/i), taking into account the varying rate of product 
formation by a factor z^(0)/p. The explicit formula is 
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As Fig. [5] illustrates, the slowing down of the reaction 
with increasing frequency at the second transition in 
Fig. [2] occurs because the high frequency of the force 
shifts the enzyme towards spending more time in the low 
efficiency conformation. The increased efffciency of the 
enzyme at intermediate frequencies sets in when the ex- 
ternal force drives the enzyme back and forth between its 
two conformations fast enough to avoid the bottleneck of 
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FIG. 6; Plot of (r) as function of driving frequency uj for 
aij,f3ij <^ kxi- The horizontal lines represent the limits de- 
rived in the text. The parameters are ku — k-u — 1, k2i — 
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The simulation comprises 2 x 10 turnovers per frequency. 



the enzyme spending long uninterrupted periods being 
forced towards the slowly reacting conformation. 

Slow switching limit. When the switching between en- 
zyme states is slow in comparison to substrate reactions, 
c^ij^Pij <IC kxi, the enzyme will mostly stay in the same 
state during a reaction cycle, and we can therefore take 
the instantaneous reaction rate to be 



1/(0) + p{2\(l))iy2, 



(9) 



where Vi = k2ikii/ {ku -\- k-u + k2i) are the fixed confor- 
mation rates and is the probability for the enzyme 
to be in conformation i at a given phase of the oscillating 
force. To find this probability we first argue in the same 
way as Eq. ([7]) was obtained to see that the (steady-state) 
probability of being in the state without bound substrate, 
given that the conformation is i, is 
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Then we note that the effective rate constant for changing 
conformation from i to j can be found as 



7y = p{E\i)a,j + {l-p{E\i))/3,. 



(11) 



This effective rate now allows the probability of being in 
conformation i to be calculated similarly to Eq. ^ as 
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_ / 7ji(0)/[7y(0) +7j-*(0)], ^ < 7y 
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where j = 3 — i. The averaged waiting time for one 
turnover can again be found as (r) = 1/9. 

Discussion. Detailed knowledge of the dynamics of sin- 
gle enzymes as well as their response to external stimulus 



is important for the understanding of biochemical pro- 
cesses occurring in living cells. We studied the response 
of a single two-state enzyme to a harmonic external driv- 
ing force in the presence of thermal noise. By combina- 
tion of analytic and simulations results, we demonstrated 
the rich response behavior of the turnover dynamics of 
the enzyme to an external force. We showed that in some 
cases there exists optimum driving frequencies that min- 
imize the turnover time of the enzyme. In these cases 
one can selectively obtain information on the reaction 
rates in the final step of the Michaelis-Menten reaction 
by choosing an oscillation in a suitable frequency range. 
Moreover, we found that the response of the enzyme ac- 
tivity can be very sensitive within a small range of phase 
(j), a signature of many biological switches. We note that 
it should be possibly to access the full spectrum of rel- 
evant driving frequencies by combining different exper- 
imental methods, such as atomic force microscopy and 
optical switching methods. 
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